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Abstract. We investigate the asymptotic properties of orthogonal polynomials for 
a class of inner products including the discrete Sobolev inner products (h, g) = 

J hg dAt + X^jLi Si=o ( c j)ff ( c i)> where [i is a certain type of complex mea- 

sure on the real line, and Cj are complex numbers in the complement of supp(/i). 
The Sobolev orthogonal polynomials are compared with the orthogonal polynomials 
^ ' corresponding to the measure (j,. 

m ' 

(N 

m . 

Q 1 ■ 1. Introduction 

On 



In the last few years orthogonal polynomials on Sobolev spaces have attracted 
considerable attention. This interest is justified for several reasons: 

(1) The comparison with the standard theory of orthogonal polynomials; 

(2) The spectral theory of ordinary differential equations; 

(3) The analysis of spectral methods in the numerical treatment of partial dif- 
ferential equations; 

(4) The search of algorithms for computing Fourier Sobolev series in terms of 
Sobolev orthogonal polynomials for norms involving the function and some 
of its derivatives. 

Despite the effort, most of the progress attained is for special types of inner 
products. A survey on the subject is provided in [2] . The results known are mostly 
connected with the formal theory: recurrence relations, location of zeros, differential 
formulas and so on. Little is known concerning asymptotic properties. Possibly, 
the only results known in this direction are contained in [15] for a special type of 
so called discrete Sobolev inner product. 
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Let fi be a finite positive Borel measure whose support, supp(/u) = S^, contains 
an infinite set of points, C M. A discrete Sobolev inner product is given by the 
expression 



(1.1) (h,g) = / hgdfi + ^J^Mj^ic^Hcj), 

J j=l i=0 

where cj e R, M jti > 0, m, Nj > (see [5]). 

By S n , n G Z+ = {1, 2, 3, . . . }, we denote the monic polynomial of least degree, 
not identically zero, such that 

(1.2) {p,S n ) = 0, peP n _i, 

where P n _i is the linear space of all polynomials of degree less than or equal to 
n — 1. It is easy to see that for every n the degree of S n is n. Furthermore let L n 
be the monic orthogonal polynomials of degree n with respect to \x. When (1.1) 
reduces to 

(1.3) (h,g) = J hgd l j, + M 1 h'(c)g'(c), Mi > 0, 

where fx belongs to the class M(0, 1) (for the definition see [19], but also (2.1) below) 
and c6l \ S^, then, as proved in [15], 

(1.4) lim = frW -<*>)', 

n^oo L n {Z) 2(f(z)(z — C) 

uniformly on every compact subset of C \ , where 

ip(z) = z + V z 2 - 1, (a/^ 2 - 1 > for z > 1). 

In [15] an equivalent expression is given for the right-hand side of (1.4). Throughout 
the paper we will use the notation <p n =4 <f, K C D when the sequence of functions 
<p n converges to cp uniformly on every compact subset K of the region D. In this 
way (1.4) becomes 

S n (z) itpjz) - p( C )) 2 

L n (z) =* Mz){z-c) ' A C(L W 

The denominators Qn of the main diagonal sequence for Pade approximants of 
Stieltjes type meromorphic functions 

j=l i=0 y 3) 

satisfy orthogonality relations similar to (1.2). In fact 

„ m Nj 

(1.6) 0= pQ n d/j, + J2J2 A ^(p( z ^(z))^l Cj , pePn-i- 
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Here Aj^ may be complex numbers, and of particular interest is the case when the 
points Cj belong to the complement of the convex hull of S^. Of course, the inner 
product associated with (1.6) is not positive definite. Therefore, a priori, it is not 
possible to guarantee that the degree of Q n is n. Nevertheless, in [8] A. A. Gonchar 
proves the following result. Assume that C [—1, 1], and \x is such that 

(1.7) *CC\[-1,1], 

L n (z) 2 



then 



A well known result (see [19]) asserts that (1.7) with [—1,1] replaced by 
is equivalent to \i G M(0, 1). In this case, consists of [—1,1] and at most 
a denumerable set of mass points contained in R \ [—1,1] which, when infinitely 
many, may accumulate only at the extreme points ±1. The best known sufficient 
condition for \i G M(0, 1) is due to E.A. Rakhmanov published in two papers [20] 
[21]. This condition is // > a.e. on [—1,1] and = [—1,1]. Simplified proofs 
may be found in [16] and [22]. 

When there is one simple pole in (1.5) (m = 1, Ni = 0) and = [—1, 1], then 
the right-hand sides of the asymptotic formulas (1.4) and (1.8) coincide although 
the corresponding inner products differ in the atomic aggregate. In the second case, 
the inner product is given by 

(1.9) (h,g) = J hgd^ + M h(c)g(c). 

The agreement is deeper than formulas (1.4) and (1.8) suggest. In [19, Lemma 
16, §7], a detailed study is carried out comparing the asymptotic behavior of the 
orthonormal polynomials (and their leading coefficients) with respect to (1.9) (with 
M > 0) with those of for all c G R. As proved in [15] these relations are the 
same if we replace (1.9) by (1.3). 

We will consider inner products for which (1.1) and (1.6) are particular cases 
(see (4.1) and (4.3)). The corresponding orthogonal polynomials satisfy relative 
asymptotic formulas of the type given by (1.8). In our approach, which follows that 
of A. A. Gonchar in [8], formula (1.7) applied at the points Cj plays a key role. Thus, 
in the following, we assume that c j G C \ and Mj^ G C and thus allow complex 
values for these quantities. We will also allow /i to be a certain type of complex 
measure which extends the class M(0, 1). These measures are introduced in Section 
3 where the asymptotic behavior of the corresponding orthogonal polynomials is 
studied. In Section 4 the main result of type (1.8) is obtained (Theorem 4). From it 
we can obtain the following result. Let S n be the nth monic orthogonal polynomial 
with respect to (1.1), where \i G M(0, 1), Mj :i G C, Cj G C \ 5^. For each fixed 
j = 1, . . . , m, we denote by N* the number of coefficients different from zero 
(i = 0, . . . , Nj). If L n is the nth monic orthogonal polynomial with respect to fx, 
then 

(L10) T^T^ll W^TT— 7) > KCC\S,. 
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Using Rouche's theorem we see that for all sufficiently large n, the degree of S n 
is n. Moreover, (1.1) implies that any neighborhood of compactly contained in 
C \ {ci, . . . , c m } has, for all sufficiently large n, exactly n — N (N = N± + ... + N m ) 
zeros of S n (for short, we say that n — N zeros of S n concentrate or accumulate on 
Sn), while each sufficiently small neighborhood of Cj, j = 1, . . . , m contains, for all 
sufficiently large n, exactly N* zeros of S n (again, for short, we say that cj attracts 
N* zeros of S n ). 

We shall limit ourselves to measures with compact support, but the methods 
employed in [10], [11] allow to extend the results in Section 4 to measures with 
unbounded support. 

Section 2 has an auxiliary character. Most of the formulas given there are known 
to the specialist, though precise reference, in the extension stated, may be hard to 
find in the literature. Therefore, for completeness, we include this short section 
with full proofs except when exact reference is available. 

2. Auxiliary results 

Let fi be a finite positive measure on KL Assume that fx G M(0, 1), that is, if 
l n (z) = T n z n + ■ ■ ■ , r n > 0, is the nth orthonormal polynomial with respect to fx 
with recurrence relation 



zl n (z) = a n+ \l n+ \{z) + b n l n {z) + a n l n -i(z), 

then 

(2.1) lim a n = \, lim b n = 0. 

n^oo Z n—>oo 

As pointed out, fi G M(0, 1) is equivalent to 

(2.2) bfiW = t (p ( z ) = z+y / z 2- li KCC\S^ 

and has the structure described in Section 1. If L n is the nth monic orthogonal 
polynomial, one also has 

(2.3) * CC \ S " 

L n {z) 2 

and 

(2.4) lim ^±1 = 2. 

From (2.2) or from (2.1) it follows that for all fixed v G Z + 

(2.5) lim l ^4^ = 2 lim = ^ z) ^ 

— l%\z) L%\z) 

uniformly on compact sets of C \ S^. See Lemma 4 in [6] where (2.5) is obtained 
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the more complicated structure of £ M for \i G M(0, 1) offers no real difficulty since 
we know that between two consecutive points in \ (— f , f), for all sufficiently 
large n, there lies exactly one zero of L n which is attracted by the point further 
away from [—1,1] (the same is true for /i^). 

In [19] (see Lemma 1.5 of Section 4.1), P. Nevai proves that for fi G M(0, 1) 

(2-6) ^^=L=, ifcC\S, 

nl n (z) y/z 2 -l 

In fact, (2.6) follows from the weaker assumption on /j that it be such that 
(2.7) IW| 1/n =^ \<f(z)\, KcC\S^. 

Obviously (2.2) and thus (2.1) implies (2.7). Also, (2.7) is satisfied if S M = [-1, 1] 
and C({// > 0}) = 1/2 where C(-) is the logarithmic capacity of the indicated set. 
This statement is contained in a little known paper of P.P. Korovkin [9]. Actually, 

(2.7) follows from a still weaker condition known as Ullman's criterion given in [26] 
(see also Chapter 4 in [25]) in terms of minimal carrier capacity or an equivalent 
criterion of Widom [28]. Orthonormal polynomials for which (2.7) holds are said 
to have regular exterior asymptotic behaviour on [—1, 1]. 

We have: 

Lemma 1. Let fi be such that (2.7) holds. Then, for all v G Z + 

(2.8) \*n\*)\ 1,n ^\<P{*)1 KdC\S„ 
and 

ll^ +1) (z) 1 



Proof. From (2.7) it follows that the zeros of l n concentrate on S^. From this and 
the fact that the zeros of l n and l' n interlace it readily follows that the family of 
functions { n-1 CAn} n>n is normal in C \ S^. Therefore, to prove (2.9) for v = 

it suffices to show that any convergent subsequence in C \ has 1 / V z 2 — 1 as its 
pointwise limit, say on M + \ S^. 
Assume that 

(2.10) lim Ij!iM=V>A(z), ACZ+, 

n^oo, n£A n L n {Z) 

for z G K C C \ Sfj,. Let z = x > Xq = max{t : t G S^}, then l n (x) > and 
<p(x) > 0. Hence from (2.7) ((2.8) for v = 0) one has 

lim ll/ n (x) = <p(x) 

n—>oo 

for x G K C (xo, oo), or equivalently 

f2.nl lim — loff l„ (x) = Ioet ud(x). 
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for x G K C (xo, oo). The derivative on the left hand member of (2.11) is equal to 
n-H'Jln. Thus, from (2.10) and (2.11) follows that 



tfjA(x) = (\oglp(x))' 



1 



and we have (2.9) for v = 0. 

Assume that (2.8)-(2.9) hold for v = k, then we show that they are also true for 
v = k + 1. From (2.9) for v = k follows, taking nth roots, that 



lim 

n— >oo 



1 li k+1) 



(x) 



n 



1/n 



= lim 



(fc+i; 



(x) 



/i fc) (x) 



1/n 



= 1, 



uniformly for a; on compact subsets of C \ S^. This, together with (2.8) for v = k, 
gives (2.8) for v = k + 1. 

Now, to prove (2.9), note that the family of functions I n _1 Z n fc+2 ' ) /^n fc+1 ' ) } is 

L J n>no 

normal in C \ S^. With this and (2.8) for v = k + 1 one obtains (2.9) for z/ = + 1 
following the same arguments as above for v = 0. □ 



Lemma 2. Let /it G M(0, 1), t/iera /or all is E Z 



(2.12) 



/ 



^n+i/(^')^n(^') 



2 — X 



dfj,(x) 



1 



Proof. Obviously, it is sufficient to consider v G Z + . By the Cauchy-Schwarz 
inequality one has for z G K C C \ S^, 



l n+l/ (x)l n (x) 



x 



dfx(x) 



< 



d(K, 5 M ) 



< oo, 



where <i(.fY, is the Euclidean distance between the two sets. Therefore, for each 
fixed v G Z + , the family of functions on the left hand side of (2.12) is normal 
and uniform convergence follows from pointwise convergence. The pointwise limit 
follows from a result by Nevai [18, Theorem 13 on p. 45]. From it we have that 



lim 

n^oo 



ln+f\X)ln\X) 



Z — X 



dv(x) = - [ 
n J-i 



1 f 1 t v (x) dx 



z-x vT 



X A 



where t u is the z/th Chebyshev orthonormal polynomial of the first kind. Thus 
(2.12) holds if we show that 



(2.13) 



1 f 1 U(x) 

TV 



dx 



1 



z-x y/l-x 2 ip u (z)Vz 2 - 1 
Recall that t (x) = 1, ti(x) = x and for v > 1 
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or what is the same 

(2.14) U +1 (x) = 2xt v {x) - t u -!(x). 

Again we proceed by induction. For v = the formula is obtained from Cauchy's 
integral formula and the residue Theorem. For v = 1 one has 

i r 1 h(x) dx z_ r 1 _j_ dx i r 1 dx 

7T 7-1 Z-X Vl - X 2 7T 7-1 Z-X y/l - X 2 ^ 7-1 y/l - X 2 

7, 

1 



1 



ip(z)Vz 2 — 1 

Now assume that (2.13) holds for v = 0, 1, . . . , k, k > 1, then we prove that it also 
holds for v = k + 1. In fact, from (2.14) and the induction hypothesis one has 

1 f 1 tk+i(x) dx 1 f 1 2xtk(x) dx 1 f 1 t}~-i(x) dx 



7T 7-i 2; -a; y/l - x 2 tt 7-i 2; - a; Vl - a; 2 W-i 

2-2 Z" 1 tfc(x) cte 2 Z" 1 t v {x)dx 



x 



2 



7r y_! 2 - a; y/l - x 2 nj-x^/l 
1 t k -i(x) dx 



X 



2 



1 r tk-i 
7T 7-i ^ - 



X y/l^X 



( 2 

2z 



V? fc " 1 (^)v / ^ 2 - 1 \¥>0*) 
1 

~ V? fc + 1 (^)v / i^T' 
which we wanted to prove. □ 

3. Relative asymptotics for certain complex measures 

As in Section 2, fx is a finite positive Borel measure, \i G M(0, 1), = supp(//). 
Let r = S/T, after canceling out common factors, where 
h t 
S(z) = l[(z - Ci ) A % T(z) = ]J(z- dj) B * , a, djEC\ A, Bj G N. 

i=0 j=l 

Set 

A = A 1 H + A h , B = Bi -\ + B e . 

Assume that Q n is the monic polynomial of least degree, not identically equal to 
zero, such that 

(3.1) = J p{x)Q n {x)r{x)dn{x), peP„_i, 

and L n the nth monic orthogonal polynomial with respect to [i. In the past few 
years, the relative asymptotic behavior of the polynomials arising from a general 
real modification of the measure \i has been considered, see e.g., [17], [18], [19] 
and [22]. We are interested in the asymptotic behavior of {Q n /L n } n G Z+. The 
corresponding result (formula (3.2) below with v = 0) was stated as Theorem 10 
in [12] but the proof was omitted being similar to that of Theorem 6 of the same 
paper. Since we also need to derive other asymptotic formulas and there are minor 
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Theorem 1. Let \i G M(0, 1). Then for all sufficiently large n the degree of Q n is 

n and for all fixed v e Z + 

(3.2) 



and 



(3.3) WU) ~ CCN 

and 

Proof. First we concentrate on (3.2) for z/ = 0. Taking p(x) = T(x)o(a;), where g 
has degree < n — B — 1, we have 



= J q(x)Q n (x)S(x)dfj,(x). 

Therefore 

A+B 

(3.5) i2n(z) := S(z)Qn(z) = A n ,fcL n+ A-jfc(2), 

fc=0 

where either A n o = 1 or Q n has degree < n. Since R n is monic, the first coefficient 
A nj fc different from zero appearing in (3.5) must be one. Dividing this relation by 

L n -B we get 

n / n Rn(z) _ v-^ x L n+A -k{z) 

n n {z :— — — > A n fe— . 

L n - B (z) ^ L n - B (z) 

Set A* = ^X^=(f l A n,fc|) < oo and introduce the polynomials 

A+B 

Pn(z) = Yl KkZ A+B ~\ P* n (z) = KPn(z). 
k=0 

We will prove that 

To this end, it suffices to show that 

/o c\ — i „ „ f„\ /A+B i \ A+B— 1 i , \ \ 
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where 



/ A+B \ _1 

(3.7) c= lim \* n = V |A fc | , A = l. 

\ fc=0 / 

Since {p%}, n G Z + , is contained in Fa+b and the sum of the moduli of the 
coefficients of p*, for each n G Z+, is equal to one, this family of polynomials is 
normal. Therefore (3.6) is obtained if we prove that for all A C Z + such that 

(3.8) lim pl(z)=pA 

n— >oo, n£A 

then pa(-z) = c Po^); where po^) an d c are defined as above. On the other 
hand, since p^ G Wa+b an d Pa ^ 0, it is uniquely determined if we find its zeros 
and leading coefficient. Finally, since the leading coefficient of p\ is positive and 
the sum of the moduli of its coefficients is one, its leading coefficient is uniquely 
determined by its zeros. Therefore p\( z ) = c Po(z) if and only if it is divisible by 
Po(z). 

Because of the factor £ in fi n and since all the zeros of L n _s concentrate on 
we immediately obtain the following A equations 



A+B , \ W 

(3.9) o = £ a;a„ J ^±^M ( Cl ), 



% = 1, . . . , h; v = 0, . . . , Ai — 1; n > uq. 



From (2.3) it follows that 



J n-B 



(i/) ( / I \ \ A+B-k\ M 



(3.10) (^2±±#) =»((^| | , KC\S, 



Relations (3.8), (3.9) and (3.10), together with the fact that </? holomorphic with 
<// 7^ in C \ [—1, 1], imply (use induction on v) 

(3.11) pM( ( rtgy\=0, i = l,...,h; v = 0,...,A i -l. 



Take p(x) = T(x)L n _ B (x)/(x - dj)", j = 1,. .. v = 1, . . . , Bj, in (3.1). Using 
(3.5) and multiplying by (z/ — l)!A*r^_ B we have the additional relations 



n 
A+B 

(3.12) = £ W- 1)!T "- B / frfr), 

j = 1,... ,£] V=l,... ,Bj. 

From (2.4) and (2.12) it follows that 
(3.13) 



T n -B 



ln+A-k{x)l n - B {x) ( -1 \(" !) 



(1/ — 1)! / 7 r (■///./• , 

V ! T n+A -kJ {X-ZY ^ } \{2^{z)) A + B - k yf^l 



is ^- m \ o 
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Relations (3.8), (3.12) and (3.13), together with the fact that l/(p is holomorphic 
with (l/<p)' ^ and 1/Vz 2 - 1 ^ in C \ [-1, 1], give (use induction) 

(3 - 14) tt'(mj) =0 ' •< : '=■• " " '* 

From (3.11) and (3.14) follows that p A (z) is divisible by p (z). Therefore (3.6)-(3.7) 
hold and 

(3.15) p n {z)=tp {z), ^CC. 

From the definitions of p n , Q n , (3.10) with v = and (3.15) we obtain 

Rn(z) ( <p(z) 



Lu-b(z) 



from which (3.2), for v = 0, immediately follows. 

We continue by induction. Assume that (3.2) is true for v = k. Note that 



r (fc+l) I r(fc) / ' ,(fc) r(fe+l)- 

Hence, from (2.9) and (3.2) for v = k we obtain (3.2) for v = k + 1. 

Using Hurwitz' theorem (see [1, p. 178]) (3.2) implies that the zeros Of Q { n\ 
for each fixed v e Z+, can only accumulate on as n — > oo. Therefore, for 
sufficiently large n, the left hand members of (3.3) and (3.4) are meaningful (finite) 
on K C C \ Sfj,. Now (3.3) and (3.4) are trivial consequences of (3.2) combined 
respectively with (2.5) and (2.9). □ 

Remark 1. If r(x) has constant sign on 1R then (3.3) and (3.4) reduce to (2.5) 
and (2.9) since then r(x) dfi(x) is a measure of constant sign on its support. So 
only the case when r has zeros or poles in a non-symmetric way with respect to M 
is of interest. 

From (2.4), (2.8) and (3.2) we can conclude that 

iQi v) (*)r /n =*^, kcc\s,. 

Theorem 2. Let n G M(0, 1). Then for all sufficiently large n (n > n ) 

(3.16) \:= [ Q 2 n (x)r(x)d»(x)>0. 

K n J 

Moreover, for such n the polynomials Q n satisfy a three-term recurrence relation 

(3.17) Qn+l(z) = {z- Pn)Qn(z) ~ a£Q n _i(z), 

where 

(3.18) 3„ = nl [ xQl(x)r(x)du(x), lim 3„ = 0, 
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and 



■2 ''" 1 i;... .-> 1 



(3.19) Q " = if' 

Finally 

(3.20) lim 4 = (-2r n '; y ' W . 



Proof. If (3.16) were not true, then for an infinite set of indices A one would have 
Qn = Qn+i for n G A, which contradicts the fact that for all large n the degree of 
Q n is n. 

To prove (3.17) we observe that the right-hand member of that relation is or- 
thogonal with respect to r(x) dfix) to all polynomials of degree < n — 2. If n is large 
enough, (3 n and o? n as given by the left hand of (3.18) and (3.19) are well defined 
(finite) and furthermore the degree of Q n -i is n — 1 the degree of Q n is n. We also 
have 



/l 
[(a; - (3 n )Q n (x) - a^Q n _i(x)] Q n _ 1 {x)r{x) d/i(x) = — ^— = 0, 
K n K n-1 

and 

J [{x - p n )Qn{x) - alQ n - 1 (x)] Q n (x)r(x) dfi(x) 

/Pn 
xQ n {x)r{x) dji{x) j = 0. 
K n 

Since the degree of (z — (3 n )Qn{z) — ^nQn-i(z) is n + 1 (the degree of Q n +\ for 
large n), we see that (3.17) follows from the above. 

To prove the limit relations in (3.18) and (3.19) we proceed as follows. Equation 
(3.17) may be written as 



zQn(z) = Q n +i(z) + f3 n Q n (z) +a n Q n - 1 (z). 
Multiplying either side by S(z) and using (3.5) we have 

A+B A+B 
(3.21) ^2 ^n,kZL n+ A-k(z) = ^2 ^n+l,kL n +A+l-k(z) 
k=0 k=0 

A+B A+B 

+ ^2 Pn^n,kL n+ A-k(z) + ^2 a n^n-l,kL n +A-l-k(z) ■ 
fc=0 k=0 

Using the three-term recurrence relation for the polynomials L n 
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the left hand of (3.21) can be written as 

A+B A+B 

(3.22) X ^zL n+A-k ,kL n +A+l-k (z) 
k=0 k=0 

A+B A+B 

+ E bn+A-k^n,kL n +A-k(z) + E a n+A-k^n,kL n +A-l-k (z) ■ 
k=0 k=0 

Since {L n }, n G Z+, is a basis, comparing the coefficients corresponding to L n+ ^ 
and L n+y i_i in the right-hand members of (3.21) and (3.22) gives 

A n +l,l + Pn^n,0 = A n i + &n+^A n n, 

(3.23) 2 _ 2 

We know from (3.15) that 

A n>0 = A n _i i0 = 1, lim A n+ i 5 i = lim X n l , lim A n+ i j2 = hm X n2 , 

n^oo n^oo n^oo n— >oo 

and since /U G M(0, 1) we also have 

lim b n = 0, lim ai = —. 

n^oo n^oo 4 

From all this and equations (3.23) we get (3.18) and (3.19). 

In order to prove (3.20), note that for n > n the degree of Q n — TL n _ B is 
< n — 1. Therefore, from this and (3.5) 

(3.24) 4 = T n l Ql(x)r(x) d»(x) = i* f L n _ B (x)Q n (x)S(x) dfi(x) 

K n J J 

2 f 2 T n 
= K,A+BT% I L 2 n _ B (x) dfJ,(x) = \ n) A+B^T~ ■ 

J T n-B 

Since 

(3.25) lim X n , A+B = X A+B =(=±} ^ (Q) 



2 J n- =1 ^K) 

(see (3.15), (3.6) and the relation preceding (3.6)), we can now use (2.4), (3.24) and 
(3.25) to obtain (3.20). This completes the proof of the theorem. □ 

Writing 

1/2 

a = ( J Ql(x)r(x) dfijx) x 
\f Q 2 n _ 1 (x)r(x)dn(x) 

(cfr. (3.18) and (3.19)), we have by (3.19) 



RELATIVE ASYMPTOTICS FOR SOBOLEV INNER PRODUCTS 



13 



Choosing uq such that 

-1/2 



we can define n n for n > no recursively by 



Kn+1 — 



With this selection we take q n = K n Qn and the sequence {q n }, n > no becomes 
orthonormal with respect to dp(x) = r(x) dp,(x), with the degree of q n equal to 
n. Moreover, because of (3.17) these polynomials satisfy the three-term recurrence 
relation 

(3.26) xq n (x) = a n+1 q n+1 (x) + f3 n q n (x) + a n q n - 1 (x), n > n , 
with 

(3.27) lim a n = -, lim (3 n = 0. 

n^oo Z n^oo 

Moreover, for all fixed v e Z+, the zeros of qn , n > no, can only accumulate on 
Sfj,, as n — > oo, and 



(3.28) 44^^^), ifcC\5 

together with 

From previous results it is easy to prove also that for all fixed v e Z 
(3.30) y |g n (^)9n+i/(aj)| |rfp(a?)| < C < oo. 

Here, as usual, \dp\ is the total variation of the complex measure p (see, e.g., [23]). 
In fact (see (3.5)) one has 

\q n (x)q n+v (x)\ \dp(x)\ = \K n K n+l/ \ [ \Q n (x)S(x)Q n+l/ (x)S(x)\- 



Hi 



S(x)T(x) 

A+BA+B dp(x) 
fc =0 m=0 7 \b{X)l [X)\ 

■ A+BA+B I \ \ I /• 

n^n+i/| \ ^ \ ^ T~n^n+j/ \ ^n,k / ^n-\-u,m \ I ; / \ 7 / 
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Set C = in£{\S(x)T(x)\ : x G S^} (> 0). Then using the Cauchy-Schwarz inequality 
one has 



A+B A+B 



I \q n (x)q n+u (x)\ \dp(x)\ < i. K n K n+v — I™ ln+Z |A n , fc A n+I , )m |, 

J ^ TnTn+iy , „ „ T~n+A—k Tn+A+v-m 

k=0 m=0 

and (3.30) immediately follows using (2.4), (3.15) and (3.20). 

In [14] A. Magnus introduced a more general class of complex measures to which, 
for large n, a sequence of (non-Hermitian) orthonormal polynomials may be asso- 
ciated satisfying (3.26) — (3.30). These are all complex measures on [—1, 1] of the 
form dp(x) = g(x)iv(x) dx, where g is a non-vanishing complex valued continuous 
function on [—1, 1] and a; is a positive almost everywhere real integrable function 
on [-1,1]. 

When f3 n 6 K, a„ > 0, the existence of a positive measure p supported in R with 
respect to which {q n } is orthonormal is guaranteed from the recursion formula 
by Favard's theorem. If, in addition, (3.27) holds (p G M(0, 1)), then p has the 
structure described in Section 1 and automatically (3.28)-(3.30) follows ((3.30) is 
trivial for positive measures), see [14], [19], [24] and [27]. 

Sequences of polynomials satisfying general recurrence relations of type (3.26) 
have been studied. See, for example, [4] for representation theorems and [7] for 
bounds on the zeros. Such relations are satisfied by the denominators of diagonal 
Pade approximants and certain types of continued fractions, hence their impor- 
tance in approximation theory. Such sequences of polynomials need not come from 
complex measures (for example, they also arise from orthogonality relations of type 
(1.6)). In case they do, what restrictions on the class of measures, where we seek 
the solution of the corresponding complex moment problem, determine a unique 
measure? This problem is commonly known as the question of determinacy. 

Recently, one of the authors in a joint paper with E. Torrano and D. Barrios [13] 
studied the location of zeros of polynomials given by (3.26) under the assumptions: 

sup(|a n |, \P n \) < C < oo, lim 3m a n = lim 3m f3 n = 0. 

n n— >oo n^oo 

In particular, from their results it follows that if (3.27) holds, then automatically 
(3.28) and (3.29) hold uniformly on compact subsets of C \ ([-1, 1]UE), where E 
is at most a denumerable set of isolated points in C \ [—1, 1] for which E' C [—1, 1]. 
Moreover, in this case (3.27) also implies uniqueness of the generating measure 
(should it exist) with support contained in [—1, 1] U E. An open problem is still 
whether under these circumstances (3.30) is satisfied. 

Definition 1. We say that p G Mc(0, 1) if the corresponding orthonormal polyno- 
mials satisfy (3.26), (3.27) and (3.30), where S p = [—1, 1] U E, and E is at most a 
denumerable set of isolated points in C\ [—1, 1]. 

Theorem 3. Let p G Mc(0, 1). Then for all f continuous on S p and fixed v G 7L 



(3.31) lim / f(x)q n (x)q n+ „(x)dp(x) = - I 

n^oo J 7T J_ x 



1 f(x)t\ v \(x) dx 



yl — ,v 



2 



where t\ v \ is the \v\th Chebyshev orthonormal polynomial of the first kind. 

Proof. If f(x) is a power of x, then the proof is carried out in the same fashion 

„„ rpi m a o i o :„ rial rpi o /o Q1 - ! l„i j„ c — „n — u. — ™;„i„ Tf :„ „ 



RELATIVE ASYMPTOTICS FOR SOBOLEV INNER PRODUCTS 



15 



continuous function on S p then, since S c p is connected and S° = 0, according to a 
theorem of M. Lavrentiev [6, Thm. 8.7 on p. 48] / may be uniformly approximated 
on S p by polynomials. This combined with (3.30) implies (3.31) for /. □ 

Corollary 1. Assume that p e Mc(0, 1). Then for all v, r\ E Z + with v > 77 



(3.32) 



n"-iqi?\z) \V^T 
and if m, fcfZ with k < m and v,r} e Z + 



v — r\ 



K <zC\S 



pi 



(3.33) 



(-1)^-1)1 /•» 



n- 



,(x)q n - k {x) 



dp(x) 
(z - x)v 



tm—k\ x ) dx 

\/ z 2 — 1 J 7T J (z — x)y/i — x 2 



1 (»»-!) 



Proof. The asymptotic behavior (3.32) follows from (3.29) and the fact that 



Let us prove (3.33) for 77 = 1. Since < m it follows that 

Qn-m( Z ) ~ Qn-m( X ) 



I 



z — X 

Therefore 



q n -k(x) dp(x) = 0. 



1_ f q y n A m {x)q n - k {x) _ q y n l m {z) f q n - k (x) 

(3.34) »"J >-* M)= ( »; ' *-* p( ' 

?n— m(^) /" ( ln—m{ x ) ( ln—k{ x 



n u q n - m (z) J z-x 

From (3.30) and (3.32) we have that, given K C C \ >Sp, this family of functions 
is uniformly bounded for z £ K. Thus (3.33) with 77 = 1 follows from pointwise 
convergence. This is guaranteed by (3.31) with f(x) = l/(z — x), (3.32) and (3.34). 
For arbitrary 77 > 1 we observe that 

1_ r q^} m (x)q n - k (x) = (-l)^" 1 ( J_ f Qn-m( x )Qn-k(x) X 

n v J "" (z-x) 1 ! " p[x) (77 -1)! In" J z-x P[X) 

and we obtain (3.33) using that formula with 77 = 1. Formula (2.13) allows to 
express (3.33) in terms of (p. □ 
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Lemma 3. Suppose that p G Mc(0, 1), rj, v G Z + . TTien 

/ 1 \ H ^ ^ 
(3-35) ^777 = 2^> fc " 



9n \ z ) / fc=0 9" l^J 

where g n ,k, k = 0, . . . , v, are rational functions whose poles accumulate on S p as 
n — > oo. Moreover 

(3.36) 9n,k=l9keH(C\S p ), KcC\S p . 

In particular, for each r\ G Z + 



^^(7^f=J ' KCC\S P . 

Proof. We proceed by induction. Fix n G Z + ; for v = the equations (3.35), (3.36) 
and the statement about g Uj o follow easily by taking g n Q = 1. Assume the lemma 
holds for v = m, and let us prove it for v = m + 1. Using the induction hypothesis 
we obtain 

! \ (m+1) / ™ n k 

,Qn"(z)J \ k=0 Qn"{z) i 

= Z^fcU) (t?)/ , - 2^9n,k- 



1 m / g^ +] 

^n,o(^)-^)— + ^ <7n,fc(*) - Sn.fc-l"^ 
9n (*) k =l V n 9n 



9n,m(z 



$\z) J Qn(z) 

qtt +1) (z) n m+1 
nq^(z) $\z) 



Using the induction hypothesis together with (3.32), the lemma readily follows for 
v = m + i. □ 

A repetition of the arguments in the proofs of Theorems 1, 2 and the conse- 
quences thereafter gives: 

Lemma 4. Let p G Mc(0, 1) and r be a rational function whose zeros and poles lie 
in C\S p . Then r dp G Mc(0, 1). Moreover, the orthogonal polynomials with respect 
to r dp and dp have relative asymptotic behavior as described by (3.2) and (3.20). 

4. Relative asymptotics for generalized 
discrete sobolev inner products 

Definition 2. Let p be a complex measure with compact support S p C C and 
ci, C2, • • • , c m G C. A generalized discrete Sobolev inner product is an expression of 
the form 

(4.1) (h,g) = / hgdp + Y.J^hVicj^iisiCj), 
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where Cj^{g;Cj) is the evaluation at Cj of the linear ordinary differential operator 
£j,i with constant coefficients acting on g and Cj,Nj ^ 0, j = 1, . . . , m. 

Since the operators in (4.1) are evaluated at a single point Cj, the assumption 
that they have constant coefficients is not a restriction, but the results hereafter 
are better described this way. The linearity is essential. 

Denote by S n , n G Z+, the monic polynomial of the least degree, such that 

(4.2) (p,S n } = 0, peP n _i. 

The existence of S n , for each n G Z + , follows from solving a linear system of n 
homogeneous equations of n + 1 unknowns. Uniqueness readily follows from being 
of the least degree. If the inner product is positive definite, then the degree of S n 
is n and thus all the S n 's are distinct. In general this is not so and different n's 
may have the same S n . 

An inner product of the form (4.1) may not be symmetric. Therefore we must dis- 
tinguish between right-orthogonal polynomials defined by (4.2) and left-orthogonal 
polynomials, which are defined by (S n ,p) = 0, p G P n -i- However, since 

n* 

(h,g}= j ghd f i + f^J2g^(c j )jCl i (h;c j ) = (g,h)*, 

J j=l i=0 

it follows that left-orthogonal polynomials with respect to (•, •) are right-orthogonal 
polynomials with respect to (•,•)*• Therefore it suffices to study right-orthogonal 
polynomials, which we will simply call orthogonal polynomials. 

Special cases of generalized discrete Sobolev inner products are (1.1) and the 
inner product associated with (1.6), that is 

(4.3) / hgd^ + J2T, h{i) ^)T, A M.)9 ik - i) (c j ). 

J j=l i=o k=i ^ ' 

Given j = 1, . . . ,m, let Jj be the maximum order of the differential operator 
Cj,i, i = 0,... , Nj. Then 

fc=0 

Let Fj = (7^ k ), i = 0, . . . , Nj, k = 0, . . . , Jj, be the matrix of the coefficients of 
the Cj^. Denote by T* the matrix obtained from Yj after deleting all the rows and 
columns with zero vectors. 

Definition 3. We say that {•,•) is a regular generalized discrete Sobolev inner 
product if for each j = 1, . . . , m the matrix V* is a square matrix with determinant 
different from zero. We denote by Ij the dimension ofT*. 

This definition means that the total number of different derivatives appearing in 
the Cj/s equals for each j the number of these operators not identically equal to 
zero, plus the fact that the non-zero operators be linearly independent. Obviously 
this is the case with (1.1) where for each j the matrix V* is a diagonal matrix, 

„„J I a o\ „.i j-U„ „ J; — , "n* „ — ,..!„„ „ „ ™„+„; — „ td„j-U „c 
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these inner products are symmetric (for (4.3) this follows easily from its equivalent 
form (1.6)), but non-symmetric regular generalized discrete Sobolev inner products 
are easy to construct. Obviously (•, •) is symmetric if for each j the matrix V* is 
symmetric. From this fact it immediately follows that not all such symmetric inner 
products are regular. 

The theorem which we are about to prove can easily be extended to certain 
non-regular inner products of type (4.1), but we restrict our attention to regular 
inner products for three reasons. Firstly, they contain the most interesting cases 
(1.1) in its complex version (Mj ;i G C,Mj :N . ^ 0, Cj G C) and (4.3). Secondly, 
the notation and technicalities rapidly increase with generality, making the reading 
(and writing) difficult. Thirdly, because we are sure that the best statement we can 
prove does not have a final character. The proof of Theorem 4 contains the main 
ingredients of our approach in solving the general problem. 

Theorem 4. Consider a regular inner product of type (4.1) such that p G Mc(0, 1) 
and ci, . . . , Cm G C \ 5^. Let {L n }, n G Z +; be the sequence of monic orthogonal 
polynomials with respect to p and {S n }, n G Z+, the monic orthogonal polynomials 
with respect to the given inner product. Then, for all sufficiently large n, the degree 
of S n is n and each point Cj attracts exactly Ij zeros of S n , while the rest of the 
zeros concentrate on S^. Also, for each fixed v G Z + 

Sn\z) TTT f {f{z) — <^(Cj)) 2 



(4.4) Tt^TT /w > KCC\S». 



Proof. Write 



>M = II(*-<*>'* +1 . ^M= ( ,_^+r 



Taking p(z) = s(z)pi(z) in (4.2), where p± is an arbitrary polynomial of degree 
< n — A — 1 with A the degree of s, we have 



(4.5) 0=( spu S n ) = J Pl S ri 



dp, 



where dp = sdp,. By Lemma 4 we know that p G Mc(0, 1), and thus, by the same 
lemma, we also know the asymptotic behavior of the monic orthogonal polynomials 
Q n , n G Z_|_, for the measure p relative to the polynomials L n , n G Z + . Therefore, 

instead of (4.4) we shall study the behavior of the ratio Sn /Qn ', n G Z+. Formula 
(4.4) for arbitrary v G Z + follows, as before, by induction on v. Thus we will restrict 
our attention to the case v = 0. Since the degree of Q n is n for all sufficiently large 
n, which in turn implies that for such n's 



/ 



Qldp^O, 



we can conclude that any polynomial S n satisfying (4.5) must be of degree at least 
n — A for n > n . Now it is easy to see that 

A 

(4.6) S n (z) = ^2\ n ,kQn-k(z), 
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since any such polynomial satisfies (4.5) and the A+l parameters X n ,k, fc = 0, . . . , A 
are sufficient to guarantee a nontrivial solution of the remaining A homogeneous 
linear relations 

A 

= {t v , KkQn-k), v = o, . . . , A - 1 , 

k=0 

which determine S n . 

As in Theorem f, the proof of the corresponding asymptotic formula for the 
sequence {S n /Q n }, n G Z + , relies on proving that for each k = 0, . . . , A, one has 
the asymptotic behavior X Uj k — > Xk and finding the limits These parameters 
are not the same as those appearing in (3.5), but in the proof they play a similar 
role. Write 

O — Sn(z) _ \^ , Qn-k(z) 



Set 



and 



Qu-a(z) ^ ' Q n -A(z) 



K= (X^lVfcl ) (<°o). 

\k=0 



hn(z) = X n,kZ A \ h* n (z) = X* n h n (z). 
fc=0 



(4.7) hn (z)^h (z):=]J(z-^) 13 (z--^) Nj+1 ^ , K C C. 



We shall prove that 

As in Section 3, this follows by proving that 

(4.8) h* n {z) =4 c /i (z) = c(z A + \ lZ A ~ l +■■■ + X A ), 

where 

hm A;= (f> fc | ) , A = l. 



c = 

n^oo 



Since the sum of the moduli of the coefficients of /i* is one, the family {/ijj}, n G Z+, 
is normal and (4.8) holds if we prove that for each subset A C Z + for which 

(4.9) hm h* n (z) = h A (z), 

the limit h\ is divisible by h . To this end, let us find a convenient system of A 
linear equations for the parameters X nj k, k = 0, . . . , A. Fix j G {1,2,... , m} and 
put p(z) = §j(z)(z - Cj) 1 , I = 0, . . . ,Nj, in (4.2). Then 

N ■ 

(4-10) 0= / S n (x) _y$ +1 _ f +Y^(l)^- I \c J )C J , l (S n ;c J ). 
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Since j is fixed we will drop it from the notation in the rest of the proof. Thus, in 
what follows we use Nj = N, Ljj = Li, Cj = c, and §j = s (which differs from the 
s used higher). Note that s(c) = §j(cj) ^ 0. This along with (4.10) allows us to 
prove that there exist polynomials pi = pjj independent of n such that the degree 
of pi is N + 1 — i, pi(0) = and 

1 



(4.11) 







= J S n (x)pi 



x — c 



dp(x) + £»(£„; c), 



i = 0, . . . , N. 



Indeed, system (4.11) is equivalent to (4.10) and is obtained from it in view of the 
triangular structure of the sum in (4.10). 

We will combine the equations in (4.11) in two different groups: those corre- 
sponding to the operators Li ^ and those where Li = 0. Let %\, . . . ,%i be those 
indices % for which Li ^ 0. According to the assumption of regularity, I = Ij is 
equal to the dimension of V* (see Definition 3). Let d±, . . . , di be the orders of the 
derivatives appearing in Li x , . . . , Li z and recall that V* is a square matrix. Since 
det(r*) ^ we can further transform the equations in (4.11) corresponding to the 
indices i±, . . . ,ii and find polynomials p\ , . . . ,p* t , independent of n, such that the 
degree of p* r is < N + 1, ^(0) = 0, r = 1, . . . , /, and 

1 X dp(x) + S^\c), r =!,...,/. 



(4.12) 







X 



In this group of equations the main role is played by the second term on the right 
hand side. In order to estimate the integral term, we proceed as follows. Take q n 
the nth orthonormal polynomial with respect to p. Since for each £ e N (the set of 
non-negative integers) 



Qu-a(c) - 



?7=0 



(c — x) 71 ) (x — c)~ 



is a polynomial (in x and c) of degree < n — A — i < n — k, it follows by orthogonality 
that 

£-1 



(4.13) 



q n -A(c) 



q n -k(x) dp(x) 
(x — cY 



D- 1 )" / 

ry=0 J 



qnl A (x)q r , 



-k 



(x) 



rj\(x — c) 1 



-v 



dp(x). 



Since the degree of p* is < AT + 1, we obtain from (4.13) and (3.33) 



(4.14) 



Qn-k(x)p* 



1 



X 



dp(x) = O 



n 



N 



qn-A(c) 



In what follows, each O(-) may be a different function of n. 

Now let us turn to the second term in (4.12). We may assume that the in- 
dices i r are taken so that < d\ < • • • < dj. Multiply both sides of (4.12) by 
\* ri K n -A/q ( ^l' A {c), where n n is the leading coefficient of the nth orthogonal polyno- 
mial q n for the measure p. Using (4.6) and (4.14) we have 

A 



= ^ Ki^n,k 



(4.15) 



fe=0 
A 



Kn-A 
^n — k 



q n -k(x) 

(d r ) / n j 



-Pi 



1 



X — C 



dp(x) + 



ai d -l(c) 



fc=o Kn ~ k 



o 



n 



N 



^q n -A(c)q^ A (c) ( 



+ 



H- ] a(c) 
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Let n G A be a set of indices for which (4.9) holds. By (4.15) for r = 1 (any 
other r gives the same result), (3.27) and (3.28) (and its consequence, the nth root 
asymptotics), one obtains by taking limits for n G A and n — > oo 



= h A 



ip(c) 



Let us now show that ip(c)/2 is a zero of h\ of multiplicity I = Ij. If the indices 
di, . . . ,di were consecutive numbers, then Leibniz' rule of derivation, via using 
(4.15), would give us this rather directly (see [8] where d\ = 0, d 2 = 1, . . . , dj = 
I — I). Here the case is a bit more complicated. To see that h' A (cp(c)/2) = we work 
out the following formula. Let p, q be arbitrary functions for which all operations 
we are about to carry out are admissible. This will be true for our functions for 
n > uq. Then 

q f p\ I V V 



q> \q) \q' q 
and thus for any r G N 

g( T_1 ) ( p^ T ~ x ^\' p( T ) p( T - 1 ) 



q(-r) \ q{r-i) J q(r) q^- 1 ) ' 
Because of the structure of the right hand side, by summing up we obtain 
p( T ) p q f p\' q' fp'\ q^ T ~^ f p^ T ~^ 



(4 ' 16) q q'{q) + q" \q' ) + "' + W T_1) 

Taking p = q^_ k , q = q^-A anc ^ t = d 2 — di, from (4.15) and (4.16) follows that 
(multiply the difference for r = 1 and r = 2 by n) 



A 

K>n—A 



= ^K^n,k- 



k=o Kn ~ k 



o 



n N+l 



^n-A{c)q^- A {c) . 



na (dl+T) (r) fa (dl+T) \ 
^ o^ 1+T+1 Vrl \ a {dl+T) / 

r=0 Qn-A \ C ) \Qn-A / 

Now, from (3.27)-(3.29) and (4.9), taking limits for n G A and n — > oo, we get 



= Vc 2 - l(d 2 - d^hA o <p/2)'(c), 

or, equivalent ly, 

V(c) 



= h' A (- 



since (fi'(c) ^ 0. This procedure can be continued until we use all the / equations 
in (4.15). We show how this is done. Leibniz' formula gives 
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Take rj from d\ up to dg + i— £. We obtain a total of di + i — £— di+1 equations. In these 
equations, the quotients p^+Z-v) / q (r)+Z-v) mn fr omp (di )j q {d 1 ) up top (d, +1 )/ ? (d, +1 )_ 

Of these quotients, we assume p^ dr ^ /q^ dr \ r = 1, ...,£ + 1, known, and the rest of 
them, unknown. In this way, we have a total of dg + i — di + l — (£ + l) = dg+i —d\—t 
unknowns. We add the fictitious 'unknown' —1, the coefficient of (p^'/g^)^. 
Put the terms corresponding to the known quotients on one side of the equations, 
then we have a system of linear equations with the same number of equations and 
unknowns. Denote by Dq the determinant of this linear system, and by _D_i the 
determinant corresponding to the 'unknown' —1, then by Cramer's rule 



(4.18) 



D = -£>_] 



Expanding the determinants Dq and _D_i along their first columns one sees that for- 
mula (4.18) expresses a certain linear combination of the expressions (p^) /q^)^\ 
77 = 1, ... , dg + i — £, as a linear combination of p( dr> ) /q^ dr \ r = 1, . . . ,£+1, say 



(4.19) 



r=l y rj=di vy 7 



Note that the coefficients j3 r and a,,, only depend on the denominator q and not 
on p (see (4.17)). Each f3 r is the sum of products of dg + i — d\ — £ + 1 factors of the 
form 



1 



7 (v) 



(tj+^-i/) 



On the other hand a v is the determinant of the minor corresponding to (p^'/q^'p^ 
in Dq. For each fixed n G A, take g = (? n -A then j3 r = (3 r (n) and a v = a v (n). By 
Lemma 3 and (3.29) 



(4.20) 



J- » "?n-A 



(?7) 



-1 



1 



9n-A,v 



l-v 



n- 



Qn-A 



+ 



(-iy 



n" 



Mn-A. 



1 



Qn-A 



rr 



KcC\S„. 



Therefore 
(4.21) 



Multiply equation (4.15) by /3 r (n), r = 1. 
formulas. From (4.19) and (4.21) we obtain 



/3 r = o(n^+ 1 " dl - £ + 1 )). 

. ,£ + 1, and sum these first £ + 1 



(4.22) = J>;A n , fc 



K n-A 
K>n—k 



o 



n 



N+e(d i+1 -d 1 -£+i) 



<ln-A(c)q^>(c) 

d e+1 -£ / ( V ) \W 
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In deriving (4.22) we recall that the coefficients (3 r {n) and a v (n) do not depend 
on the numerator p, therefore (4.19) is also valid with the same coefficients for 
p = q n -k, k = 0, . . . , A. From (4.20) we also know that a v (n), for each r/, is of the 
order n i ( de+1 ~ dl ~ i ' . Therefore, dividing (4.22) by this quantity, taking limits for 
n -> 00, n E A, we obtain from (3.27)-(3.29), (4.9) and (4.20) 



= 



lim 

n— >oo, n£A 



di+i — 

E 

V=di 



«(d £+1 -di-£) 



n 



(*) / ¥>(c)\ _ 



A 



(Vc 2 - 1)^+1 - d i-^) ' 



where M 7^ is the determinant of the matrix obtained by deleting columns 1, 
<&i — d\ + 1, . . . , d^+i — d\ + 1 (counted from right to left) of the matrix 



di + i — d\ — I + 1 rows 



/I 
1 








k — 1 columns 

-(D 
(S -0 



(-D<G)\ 


y 



M Q ... ... 

(see (4.17)). This matrix is the result of shifting the row ( ), — (J, . . . , (— to 
the left k — 1 times and padding the remaining part with zeros, but with a column 
of one's on the extreme left. Since this determinant is nonzero it follows that h\(z) 
is divisible by (z — (p(cj)/2) I i for each j = 1, . . . , m. 

Now we use the N + 1 — / equations in (4.11) for those indices i for which £j = 0. 
Before doing this let us look more closely at equations (4.10) in order to derive one 
more property of the polynomials pi. If one writes down the equations (4.10) one 
finds (again we drop the index j) 



= S, 



dp 



x 



- + M-£ N 
c 0! 



0= S 



dp 



= S, 



I 



+ m^C + (N- ly.^C^ + (N- 2)!^£ K _ 2 
do s" s' 



(N-3)\-L N . 3 



Sr. 



(X - C)^ 1 



S (JV) 



(iV^l)! 



£iV-l + 



+ l ! £ I £ 1+ 0!^£„. 



From this, finding £>tv is immediate. To obtain pn-i we take (3[ = 1 and /? 2 
satisfying 
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and from the first two equations obviously 



0= / s. 



p'l 



(x — c) ; 



+ 



x — c 



dp+(N -l)\sC N - U 



from which pn-i is determined. Taking P[ and (3' 2 as above, we then take f3' 3 such 
that 

and from the first three equations one sees that 



= S, 



P'l 



+ 



ft 



(x — c) 3 (x — c) 



from which we can find Pn-2- Continuing this process we obtain all pi (and 
■ ■ • 7 P'n+i)- 

Let i\ < ■ ■ ■ < tN+i-i denote those indices i for which d = 0. Note that 

degree p ir = d' r = N + 1 - i r > 1. 

Therefore all p ir are of different degree and their coefficients do not depend on n 
and moreover their general form is such that 



(4.23) 



Pi 



1 



x — c 



E 



As before, we show how the proof works in the first two steps, and then give the 
general outline on how to carry on. If we multiply (4.11) for i r by q^ A dr \c) and 
use orthogonality (in the second equality we use the same arguments as those we 
used to obtain (4.13), and then in the third equality we eliminate other terms also 
equal to zero by orthogonality) we obtain for r = 1, . . . , N + 1 — I 



r 1 I dpix) 

X — C 



(4.24) _ / Sn ^ E 

J v =o 

8=1 J 



<-l (d{-d' r +T,) 

1n-A 



7] 



(c - x) v p ir 



x — c 



d'r-S id'i-d'r+v) 
n^n-A 



Em)' 

?7=0 



dp(x) 
1 



T]\ (X - c) d 'r + 1 - S - r > 



dp(x). 



Now it is easy to prove that if A C Z + is such that (4.9) holds, then 

1 



0. 



In fact, note that in (4.24) the highest order of derivative appearing is d[ — 1 (for 
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by \* n K n -A/n d ' 1 1 , using (4.6), (3.33) and (2.13), we obtain by taking the limit for 
n — > oo, n e A, that 



1 



4-1 



h A 



1 



we need to eliminate the term with derivative 



In order to prove that h' A y 2 (p(c) 
d'i — 1. In order to do this we combine the first two equations. Note that d[ > d' 2 > 1 
(actually d 2 > 2 because Cn ^ 0). Since 



K-l)! 

1 



(4-2)! 
(4=T]T (4=2)1 



there exist a± and «2 such that 

(_l)<-i 



and 



+ a 2 



(-1) 



4-1 



K-l)! '(d^-l)! 
+ a 2 - 



= 1. 



K-2)! ' ^(4-2)1 

Multiply the first equation in (4.24) (r = 1) by ai\* n K n - a / n d>1 ~ 2 , the second one 
(r = 2) by CK2A* K n _^/n dl_2 , sum them up and take the limit for n — * 00, n e A. 
Again using (4.6), (3.33) and (2.13), we obtain (even if (3 2 7^ because h\ y 2ip ( c ) 
0) 



= 



1 



1 



4-1 



2p(c) 



Using h\ ^2^c)) = ^ this f° rmma implies 



2^(c) 



The proof is completed by induction using the scheme employed in passing to the 
first derivative. Of course the situation becomes more complicated since each time 
there appear more terms with the same order of derivative. Fortunately, this causes 
no real problem because the polynomials pi r have correspondingly equal coefficients 
starting from the highest degree down (see (4.23)). To be more precise, in proving 
that 



h 



(t-1) 
A 



1 







T<N+1-I, 



we start with the induction hypothesis that 
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Noting that d[ > d' 2 > ■ ■ ■ > d' T > 1 (actually > 2) we will show that 

li l 



K- 
i 


1)! 


K-2)! 


K- 


1)! 


(4-2)! 


i 




1 



K-r)! 



K-l)! K-2)! (d' T -r)\ 

Indeed, this statement is obviously equivalent to 



1 d[-l {d[ - l){d[ - 2) 
1 4-1 (4-l)(4-2) 



1 <-l «-!)«- 2) 



This determinant can be transformed into a Vandermonde determinant by elemen- 
tary column operations. Therefore, there exist at, . . . , a T such that 



r=l 



("I) 



d'-k 



{d' r -k)\ 



and 



= k = 1,... ,r- 1, 



^ '«-t)! 



r=l 



Take the first r equations in (4.24), multiply the rth equation by ct r A* « n _^/n dl_T , 
sum them up and observe that in the sum all derivatives of order d^ — 1,... , d[— r+1 
cancel out. Using (4.6), (3.33) and (2.13) and taking the limit for n — > oo, ti£ A 
then gives 



o = £# 



r=l 



1 



<-T+l 



1 

2^ 



(r-r) 



(c), 



and using the induction hypothesis we immediately obtain that 



h 



(r-l) 



1 



N+l-I 



0. 



Hence /ia(z) is also divisible by iz — 2 ^{ c ) ) ^ or eacn J = 1? 

(4.7) holds and thus 



SJz) 



n 



<p(z) <p(cj)\ Ij ( <p(z) 1 



Q n - A (z) fj; V 2 2 7V 2 2^(cj) 

or, equivalently on account of (3.3) (see Lemma 4) 



Ni + l-Ii 



m. Therefore 



KcC\S„, 



ud(z) 



1 - 



1 
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Finally, from Lemma 4 we also have that 

It is easy to verify that 

'ip(z) - ip(c)\ / 1 



2(z-c) / V ¥>(*Mc) 
hence (4.25) and (4.26) give 



= 1, z,ceC\[-l,l] 



which is the same as (4.4) (for v = 0). For arbitrary z/ as in Theorem 1, (4.4) follows 
by induction. The statements concerning the degree of S n and the asymptotic 
behavior of its zeros follow from (4.4) and Hurwitz' theorem. □ 

Formula (4.4) and the existence of the limits of the coefficients X nj k allow us to 
obtain, using arguments similar to those given above, other asymptotic formulas 
for the polynomials S n . We collect some of them in the following corollary, the 
proof of which we leave to the reader. 

Corollary 2. Under the hypothesis of Theorem 4> f or oil sufficiently large n, 
(S n , S n ) >t^ 0. Let r n be the leading coefficient of l n , the nth orthonormal polyno- 
mial with respect to /i. Then 7 n = (S n , Sn) -1 ' 2 , n > no, may be taken so that 

m, 

In 



lim — 1 | 



and in particular 

lim ^±1 = 2. 

Denote by s n = 7 n Sn ; n > no, the nth orthonormal polynomial with respect to (■, ■). 
Then for all fixed v e Z + 



s 



ls^ +1 \z)_ 1 



KcC\(^U{c 1; ... ,c m }). 



Theorem 4 and the asymptotic properties of orthogonal polynomials with respect 
to measures in M<c(0, 1) allow us to obtain an extension of Gonchar's result in [8] 



2<S 
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Corollary 3. Let f be as in (1.5) where fi E M<c(0, 1). Let n n be the [n — l,n] 
Pade approximant for f , Q n its denominator with leading coefficient equal to 1 and 
L n the nth monic orthogonal polynomial with respect to ft. Then 

For all sufficiently large n the degree of Q n is n, and each cj attracts exactly Nj + 1 
zeros of Q n , while the rest of the zeros of Q n accumulate on S^. The rate of 
convergence of n n to f is given by 



f(z) - n n+1 (z) _ 1 



f(z)-n n (z) p 2 (z) 

In particular 



KcC\(S M U{ci,... ,c m }). 



lim ||/-7r n ||£ = -n — □ — , ifcc\(5 (1 u{c 1) ... , Cm ». 

n^oo \\^P\\K 

As a final remark we wish to say that for regular generalized Sobolev inner 
products it is not difficult to prove that the zeros attracted by the c/s converge 
geometrically when /j, 6 Mc(0, 1). Special cases show that in general this is not so 
if the inner product is not regular. Nevertheless, it is plausible that formula (4.4) 
remains valid in the nonregular case if we substitute Ij by the rank of V* . 
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